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We study the line shapes of radiative (/!>-decays involving virtual /o(980) and ao(980) mesons which 
decay, via derivative couplings, to tt^tt" and respectively. After developing the formalism for 
derivative interactions at one-loop level, we show that they can reproduce the measured peaked line 
shapes of 0-decays without including kaon loops. 



I. INTRODUCTION 

'66'. 

■ The radiative decays of the (f> meson are a valuable tool to study the nature of light scalar states below 1 GeV 
! P) S I3l ■ It is still deljated if the latter are quarkonia [1] , tetraquark B , molecular states or a mixing of these 
configurations [1, Q . A dominant quarkonium assignment is problematic , thus leading to identify the gg-states 
' with resonances above 1 GeV where mixing with the scalar glueball takes place [ll|, [l^l • 

. In this article we do not employ a particular interpretation for the resonances /o = /o(980) and = ao(980). 
Independently on their nature, spontaneous breaking of chiral symmetry implies that derivative couplings to the 
\l , pseudoscalar mesons tt, rj and K arise. In fact, the latter states are the emerging Goldstone bosons, which can 
' ' be rewritten as angular excitations. Then, Lagrangian interactions of derivative type with schematic form Lint ~ 
,— S {d^ipi) {d^ip2) with S — /o, ao and ^pi^2 — tt, rj, K, are obtained [l^l- Our aim is to study in detail the effects of such 
' derivative interactions on the line shapes of the reactions 4> ~* 77r''7r'' and (j) ~* Jt^^Vi occurring via virtual /o and oq 
Q-i' mesons respectively. 

^ • In order to discuss properly the issue we concentrate on the reaction (j) — > jtt^tt'^ via a virtual fo state. We first 
^ ' write down explicitly the general interaction Lagrangian of the fo meson with pions and kaons as emerging upon 
, chiral symmetry breaking within the nonlinear realization of chiral symmetry as obtained in Ref. [isj : 

m ■ ^»t./o = c/„,,/o {df^lff + d/„,,M2/o¥^2 ^ Cf.KKfo {d^K+) (d^K-) + df,KKMlh{K+K-) + ... (1) 

^sj where dots refer to the analogous terms with the neutral kaon states. In the chiral limit the pion and the kaon 
masses A^^^ and vanish leaving only interactions with derivatives parametrized by the coupling constants Cf^-^-^ 
and CfgKK- The constants d/oTTTr and dj^KK parametrize the interaction without derivatives. Now, the radiative decay 

" [ , (j) — > 77r''7r° via fo can occur essentially in two ways: 

(a) An interaction Lagrangian for the decay (j) 7/0 is introduced as C^jf^ = c^^/o/o-F^'^V^i/, where F'"' = 
, d^A" — d^Af^^ is the electromagnetic field strength and V^,y — df^tpv ~ d^tjin is the field strength related to the 

vector field 0^ . The interaction Lagrangian C^-yf^ , parameterized by c^-yf„ , corresponds to a point-like coupling which 
effectively takes into account loops of quarks (whose precise form depends on the microscopic interpretation of the fo 
meson). After the transition —^ 7/0 the scalar meson fo decays into pions (as depicted in Fig. l.c) via derivative 
couplings (mechanism a.l, whose amplitude is proportional to c/qtttt) or via non-derivative couplings (mechanism a. 2, 
proportional to df^T^-n). A first study of derivative interactions (thus setting d/oTTTr = df^KK = 0) has been performed 
in Ref. where it has been shown by using Breit-Wigner propagators for the scalars that derivative couplings 

describe the data better than non-derivative ones. A detailed study involving non-derivative interactions only (setting 
c/oTTTT = CfoKK = 0) has been performed in Ref. [Tsj . 

(b) The vector meson </> couples strongly to kaons. Then, via a kaon-loop a photon is generated at one vertex and 
a /o meson at the other vertex. The latter coupling can again occur in two ways: via derivative coupling (mechanism 
b.l, proportional to Cf^xK, see Ref. [3]) or via non derivative coupling (mechanism b.2, proportional to djgKK, see 
Ref. @,[lllI3). 

Note that all these considerations apply also for the ao meson as a virtual state of the reaction (j) ^tt^t]: the 
derivative (non-derivative) interactions with nr] and KK are parameterized by CagTiT] {dao-Kri) and CagKK (daoKK) in 
a Lagrangian which is analogous to Eq. ([T]). Clearly, the possible decay mechanisms are also separated into a.l and 
a.2 (direct, non-structure couphng) and via kaon loops (b.l and b.2). 

It is evident that the description of the radiative decay (p 77r°7r° (and (j) jir'^i]) is difficult because four decay 
mechanisms (a.l, a.2 and b.l, b.2) can potentially contribute and it is not clear a priori which one is dominant. As 
a consequence, considering that the inclusion of all contributions at the same time has not yet been performed, the 
extraction of the parameters of Eq. ([T]) from experiments depends also on theoretical assumptions. 

Mechanism (b) is a mesonic 1-loop contribution to the decay mechanism, which is regarded as dominant by many 
authors in view of the large coupling to kaons of the vector meson </> meson and of the scalars ao(980) and /o(980). 
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However, notice that the mechanism (a) is dominant according to large-Nc counting rule both in the tetraquark and 
quarkonium assignments for the /o(980) meson and in the tetraquark assignment for the ao(980) [ll] (see also the 
note to this Reference). While it is not yet clear to which extent large-Nc is reliable in this context -see also Ref. Q 
for a discussion of this point-, we consider this fact as a motivation to study in detail the effects of the mechanism 
(a) on radiative decays and thus to test an alternative scenario for the description of radiative (j) decays. Moreover, in 
this work we restrict to the chiral limit dominant mechanism a.l where only derivatives are involved (correspon ding 
to setting df^T^-^ = df^KK = daa-Kri — daaKK = 0). We thus intend to continue the analysis initiated in Ref. |14l | 
about derivative interactions by studying the decay mechanism a.l in both the /o and the ao channels in relation to 
radiative </) decays. We aim to do it by properly taking into account loops and finite-width effects using the formalism 
developed in Ref. [HI extended to the case of derivative couplings. That is, both real and imaginary parts of self- 
energy contributions (which show a rather different behavior than their non-derivative counterparts, see Fig. 2) are 
taken into account. We regularize the model by using an effective cutoff of the order of 1 GeV which is introduced by 
using a nonlocal extension of Eq. ([T]) . We conclude this discussion by noting that the neglect of the non-derivative 
coupling is surely justified in virtue of the small pion mass in the tttt channel of /o but is less justified in ttt; channel 
of ao and in the kaon-kaon channel of both resonances because of the larger masses of the K and 77 mesons. The next 
step shall be the inclusions of mechanism a. 2. but at the present stage a fit with non-zero c's and d's at the same 
time would not be constrained enough. 

Quite remarkably, our analysis shows that derivative interactions alone work well in the description of both cj) — > 
77r°7r° and '-fn^r] line shapes as experimentally measured in the SND and KLOE collaborations in Rcfs. 0, Hlj 
for the /o meson and in Refs. [13, for the ao meson. The peaked line shapes can be reproduced in virtue of the 
derivative coupling which enhance the theoretical curves close to threshold. 

We also compare our results with Ref. where the /o meson has been studied studied in j/^ decays at BES, 
and with Ref. 25], in which an experimental analysis of the Crystal Barrel data for the ao meson has been performed, 
see also Refs. |26i, i27|]. However, we shall also point out that care is needed when extracting the coupling constants of 
scalar-to-pseudoscalar from radiative cj) decays alone because of strong correlations of the parameters entering in the 
fit. In particular, the difficulty is due to the the fact that, if the decay mechanism (a) is dominant, an extraction of 
CfgKK (and eventually df^Kx) from the line shapes of — > 77r''7r'^ reaction is hard because Cf^xK and df„KK appear 
only in the denominator of the propagator of /o and are not directly proportional to decay amplitude(s). (This is not 
the case for mechanism b, where the kaon-loop amplitudes are directly proportional to Cf^xK and df^KK)- A similar 
discussion holds in the ao case. For this reason we did not include the coupling to kaons as free parameter of the fit 
but we repeated the analysis for different values of the latter checking how the fit is affected. 

As a comparison we also perform the fit using non-derivative couplings only (that is we set c/gjrTr = CfoKK = and 
we leave d/oirTr and dfgKK free) and we show that a bad description of data for the fo meson is obtained. In the ao 
case {caaTrr] = CaoKK = and daairr) and da„KK frcc), whilc a fit without derivatives is still acceptable, the resulting 
parameters are not compatible with the experimental results of Ref. [25j . Thus, the need of including derivative-like 
couplings is favoured by our study. 

The paper is organized as follows: in the next Section we extend the formalism of [19] to the derivative interaction 
and we derive the theoretical expressions for the lines shapes of 0-decays. In the third Section we present the fits to the 
experimental results and we study the correlation of the parameters. In the fourth Section we drive our conclusions 
and outlook. 



II. (t) DECAYS WITHIN DERIVATIVE INTERACTIONS 
A. 1-Loop within derivative interactions 

We generalize the study of Ref. 19] to which we refer for a careful treatment of the definitions and relative 
discussions, by considering the following nonlocal interaction Lagrangian of derivative type: 

^ = \{dt^Sf - ]^MlS^ + \{d^^f - imV' + A„t ; A„t = 9S{x) j A^dM^ + V^d^^ix - y/2)$(y) (2) 

with = d/dx'^ . The nonlocality describes the finite dimensions of the scalar states [l^: it takes into account already 
at the Lagrangian level of a form factor in the expression of the decay width. Previous studies (see, for instance, 
Refs. [Hl,!!^ and also the more microscopic approach of Ref. [l^) show that a cut-off of the order of 1 GeV emerges 
in the context of phenomenological mesonic theories. By introducing the function /a (9) as the Fourier transform of 
$(2/), ff^{q) = / (i'*y$(y)e~*^*, the tree-level Feynman amplitude of Fig. [T]a reads —ig [qi ■ 92) /a ((91 — 92)/2) where 
qi and 52 are the momenta of the two particles Lp. Let p = gi + 52 be the momentum of the particle S. 
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(a) Nonlocal derivative vertex (b) Self-energy (c) Radiative tp-decay 



FIG. 1: Relevant Feynman diagrams [colors online]. 



The decay width is evaluated in the reference frame of S, in which we have: 



(3) 



The tree-level decay reads explicitly: 



r 



{V2g) 



M2 - 2m2 



n 2 



/a(0, 9^) 



e{Mo - 2m) 



(4) 



where the equality qi ■ q2 — (Mq — 2m^)/2 has been used. The tree-level propagator As{p) = [p^ — A/q + ie] ^ , 
valid in the limit 5 — > 0, is modified by the 1-loop correction and takes the form A.s{x — ~ 
[x^ - M2 -I- (\/25)2e(x, to) + ie\ in which the self-energy m), see Fig. [TJb, is given by: 



'p^, to) 



A [(gl-g2)/A(g°,^)]' 

(27r)'* [qf — + ie] [gl — + ie] 



(5) 



and of 



with qi = q+ p/2 and (72 = pI2. For future convenience we consider the loop as a function of a; = • 
the mass to. Note that the only difference with respect to the non-derivative study of Ref. [19] is the extra-factor 
{r\ ' I"!) in ths tree- level decay rate and in the numerator of the self energy. 

A possible choice of a three-dimensional cutoff, which makes the model finite, corresponds to ( [HI, [H, [2^ ) /a((?) = 
/aCq*^) = cxp[— "g'^/A^] with qiptpix) — \/~q^ as in Eq. ([3]). The precise form of the cutoff function does not influence 
the results as long as convergence is achieved. Moreover, despite the fact that the superficial degree of divergence of 
the integral in Eq. (O in the limit of large A is four, it turns out that the EqM^ is only linear divergent in A. In our 
work we use A = 1.5 GeV, which is very close to the value used in Refs. 0,0]. Variations of the latter between 1-2 
GeV affects only slightly the results. 

A general property for E(j3^) follows from the optical theorem: 



Is(x) = (\/2g)'lm[I](a;,m)] = xV%^{x). 



(6) 



The imaginary part of the self-energy diagram is zero for < a: < 2m and nonzero starting at threshold. The real 

part Rs{x) = (V2.g) Re[S(a;,TO,)] is nonzero below and above threshold. In Fig. [2] the functions Re[S(a;, to)] and 
Im[I](x,TO)] are plotted as an illustrative example. 

We define the (nominal) mass Ms for the scalar field S as the solution of the equation M| — AIq + Rs{Ms) = 0. 
When the function R{Ms) is positive, which is usually the physical case (Fig. 2), the dressed mass Ms is smaller 
than the bare mass Mq, showing that the quantum fluctuations tend to lower it. We consider the case Ms > 2m, 
thus no pole below threshold is found and we have a truly resonant state. The spectral function ds{x) of the scalar 
field S related to the imaginary part of the propagator is: 



ds{x) = — 

IT 



lim Im[A5(a;) 

e— >0 



Ms>2m 2x 



Isix) 



IT (x^ 



M^ 



Rsix)r + isixy 



(7) 
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FIG. 2: Real and imaginary parts of the self-energy E(a;) for the illustrative values m = 0.5 GeV and A = 1.5 GeV. 



In the limit g — )■ wc obtain the desired spectral function dsix) — S{x — Mq). The normalization of ds{x), i.e. the 
validity of the Kallen-Lehman representation 30], holds for a large range of g: /p ds{x)dx — l.The decay rate for 
the process S i^ip, which includes finite-width effects, can be defined as Tsipip = /q°° dxds{x)T^g\p^{x) . This formula 
reduces to the tree-level amplitude rgJ,^(Mo) of Eq. (jH) in the hmit of small g: rgJ^^(Mo) ~ Tsip^p for g — > 0. We shall 
however not use this formula in the present work: in fact, while mathematically correct, it presents some practical 
mismatch in the case of derivative couplings because a long-tail of the function ds{x)T'''g'^^{x) may arise at large x. 
We shall therefore not compare integrated decay widths but decay amplitudes evaluated on shell, which are free form 
these ambiguities. 

As a last step of this subsection we study the case of derivative interactions of the scalar field S with two different 
particles Lpi and Lp2 with masses mi and m2'. Cmt = gSd^(fid^Lp2 and its nonlocal extension. The previous formulas 
change as follows: 



M§ - ml- m4 



/A(0,t) 



9{Mo — mi — TO2) 



(8) 



where gi^upjl-^^o) — 2M^\/^o + ("^1 ~ ^2)^ — 2(m^ + m2)MQ. The corresponding loop function is denoted as 
S(x,r7ii,TO2) and reads 



T,{x, mi, m2) 



dU 



[to • 92) /a (9°, 9^)]' 



(27r)4 [^2 -ml + ie] [q^ -ml + ie] 



(9) 



B. The /o(980) case 

1. Derivative interaction with pions and kaons 

The field /o, describing the resonance /o(980), interacts via derivative couplings with pions and kaons. In the local 
limit one has: 



C.ntjo = {d^^f + Cf„KKfo {{d^K+) (d^'K-) + {d^K^) (a^lf)) 



(10) 



where c/^tttt and Cf^KK are the coupling constants. The nonlocal case, which is used in the following calculations, is 
obtained by delocalizing the previous interaction Lagrangian: 



C^nt.Jo^gh[x) d^ydf,lr{x + y/2)d^'lr{x-y/2)^{y)+ 'kaon part'. 



(11) 
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The tree-level decay formulas, as function of the mass x, read: 



foKK 



(x) = ^f^f^ [AfoKKix)f e{x - 2mK); Af„KKix) = V2Cf„KK 



(12) 
(13) 



where for future use we introduced the amplitudes Af gT^T^ix) and Af^KKix). The field /q is dressed by pions and 
kaons (Fig. [T]b). Thus, the real and imaginary terms of the self-energy include contributions from both loops: 



Ifoix) (VQcfg^rr^ Im[i;(x, m^)] -I- (^V2cfoKK^ Im[S(a;,mx)], 
^fo{x) ^ (VQcfoyrTT^ Re[j:{x,m„)] + (y2cfoKK^ Re[T,{x,mK)]. 



The optical theorem holds for the single channels 

2 



^ lm[T,{x,mT,)] ^ xT^f~^^^{x), (^V^Cf^KK^ ^T^i^ix^rriK)] = xT^f^i^j^{x). 
The spectral function df„{x) is just as in Eq. ([7]) upon setting S — /o- 



(14) 
(15) 

(16) 



2. — > 77r'^7r'' decay 



Be F'^'^ = d^A'^ — d'^ A^^ the electromagnetic field strength and V^^, = d^(l)i, — the field strength related to the 
vector field 0^, which describes the resonance (/)(1024) of the PDG tSlJ. The Lagrangian which describes the process 
(j) — + '-yTT^TT*' and the corresponding tree-level decay rate read 



(ml - x^) 
247rm?. 



(17) 



When /o is on shell one sets x — Mf^, However, we are interested to the subsequent decay of /o into Tr^Tr*^, in which 
/o is a virtual state as depicted in Fig. [T]c. The partial decay rate, defining the line shape of the 77r°7r° decay, 
reads 



dr^jT^o^o (x) 



dx 



07/o 



(x) 



2x 



xT 



[x^-Ml^^+Rj,ix)) +If,ixy 



(18) 



where r^^'^o^o(x) — ^r^^^^^{x). It can be also rewritten as 

dr^^^o^o{x) = r*^-^^^(a;) [df^{x)dx 



^foTvO^roix) 



r*;'_(x)+r*;VK(^) 



(19) 



whose interpretation is straightforward: r^^y^(x) describes the decay rate for (j) — s- 7/0, dfg{x)dx represents the 
probability that the particle /o has a mass between x and x + dx and finally T^j^^^oT^oix) / {T^f^^„{x) +rjj^^(a;)) 

describes the branching ratio of /o decaying into a tt'^tt'^ pair. The quantity dF^^^o^o (a;)/cix can be directly compared 
with experiments as we do in the next section (32| . 



C. The ao(980) case 

1. Derivative interaction with ix^rj and kaons 

The discussion concerning the resonance ao(980) follows the same line of the previous subsection. The interaction 
Lagrangian for the neutral field Oq = aQ(980) reads in the local limit 



= '^ao.A {d,V) + Ca^Kxal ({d^K+) {d^K-) + {d,K") {df^K") 



(20) 
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where the coupling constants CagTrri and CagKK have been introduced. The tree-level decay rates in the nonlocal case 



are 



r*;W(^) = [AaoKKixf e{x - 2™^); ^..^^(x) = V2c,„kk ^ /Afe/.(x)) (22) 

where the amplitudes AaQTrrjix) and A^oifif (x) have been introduced. The real and imaginary parts of the loop include 
contributions from 7r?7 and KK loops: 

Iao{x) = (cQo^^)^Im[S(2;,TO^,m^)] + (^V2caf,KK^ Ini[i;(a;, m^)]; (23) 
■Rao(a;) = (cqottt,)^ R.e[S(x,TO7r,TO^)] + (^\/2cQoK/f j Rcp (x, m_R-)] . (24) 
The spectral function daa [x) reads as in Eq. ([7|) upon setting S — a^. 

2. 4> 77r''r; decay 

The interaction Lagrangian describing the process (f> — > •yir'^rj and the corresponding tree-level decay rate are given 
by: 



(ml - x^) 
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where x = Ma^ for an on-shell decay. However, the Oq meson decays subsequently into ir'^r] as depicted in Fig. [T]c. 
As a result the line shape of the reaction (f> — > jn'^jj reads: 



dx -'^TaoW 



[~ (:r2-Af„%^_+i?„„(a;))' + /,,(x)2 



(26) 



III. FIT OF THE LINE SHAPES 

A. The 4) ^ 77rV° fit 

Five parameters A, Mfg,c^jfg,CfgTrTr,CfgKK determine the (f> jir'^n^ line shape via a virtual /o(980) meson [33| . 
We fix the cutoff A = 1.5 GeV; a mild dependence of the results is seen by varying A. Being the coupling to kaons of 
particular theoretical interest we perform the fit for different values of cj^kk- We start setting CfgKK = ; the fit is 
done to the experimental data of SND and KLOE collaborations (20l[2l| corresponding to the black and grey dots in 
Fig. [3] retaining the data points above 0.6 GeV. The solid line is the corresponding theoretical curve for cj^kk — 0, 
whose fit parameters are listed in the first entry of Table HI The peak of the /o line shape is a feature which is easily 
reproduced by derivative interactions. 

Let us now investigate how the fit changes by increasing Cf^xK from to a maximum value of 14 GeV""'^. As we 
can see in TablelH the fitted mass varies slowly from Mf^ ~ 981 MeV to Mf^ ~ 971 MeV, c/qtttt increases from 1.28 to 

2.90 GeV~^and the also increases from 1.7 to 4. We also present the on-shell amplitudes Af^T^j^ = ^/^^^(M/q) 

and AfgKK = AfgKKiMfg) defined in Eq. (fT2| -(|13 p . which in the following will be compared with the results of [2^ . 
The value of /d.o.f increases by increasing Cf^xK- However, as already mentioned in the Introduction, the coupling 
constant Cf^^KK appears only in the denominator of the propagator of the /o meson (see Eq. (jlSp ) and therefore 
its determination is difficult. Nevertheless, we notice that rather acceptable fits are found when varying Cf^xK in 
such large range. We thus not intend to state that c/^kk is small (what would be in conflict with data from tttt 
scattering) but only that its determination from the line shape is problematic when mechanism a. 1 (whose amplitude 
is proportional to Cf^T^j^) is regarded as dominant. 
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FIG. 3: Branching ratio "^^d"/ [MeV"^] as function of the invariant mass X [colors online]. — 4.26 MeV is the 

full width of the (j) meson. We consider data sets from the SND and KLOE collaborations [20l . [2H corresponding respectively 
to the black and grey (green online) dots. The continuous line is the result of the fit by setting C[^kk = 0, the dashed line 
corresponds to the case CfgKK = 12 GeV~^. Both cases are in Table I. The dotted line corresponds also to c/^kk = 12 GeV~^ 
but only data points above 0.8 GeV are used in the fit, see Table II. 



TABLE I: Fit of the (j) ^ ^ir°-K° line shape data from x > 600 MeV. 



CfoKK 

(GeV-') 


(MeV) 


C/OTTTT 


Ctl>-ffo 

{GeV-'') 


d.o.f. 


AfgKK 

(GeV) 


A-fg-n-K 

{GeV) 





981.8 ±2.8 


1.28 ±0.04 


0.249 ± 0.008 


1.7 





1.31 ±0.042 


2 


981.7 ±2.8 


1.30 ± 0.04 


0.260 ± 0.008 


1.9 


0.675 ±0.01 


1.34 ± 0.045 


4 


981.7 ±2.8 


1.38 ±0.05 


0.291 ±0.010 


2.3 


1.35 ±0.01 


1.42 ±0.06 


8 


974.4 ± 2.6 


2.00 ±0.06 


0.312 ±0.04 


3.4 


2.62 ±0.03 


2.02 ± 0.06 


12 


971.8 ±2.3 


2.59 ± 0.06 


0.383 ±0.011 


3.9 


3.89 ± 0.04 


2.61 ±0.06 


14 


971.2 ± 2.3 


2.90 ± 0.07 


0.425 ±0.11 


4.0 


4.53 ±0.04 


2.92 ±0.07 



Notice also that at low x the resonance /o(600) contributes to the total branching ratio as found in Refs. [17ll26|. 
For this reason we repeated the analysis by retaining only the data points above 800 MeV where /o(600) is less 
relevant. The corresponding results are presented in Table [III where significantly smaller are obtained. Also 
for large values of Cf^KK the results are satisfying. As we increase CfgKK from to 14 GeV~^, the fitted mass is 

almost unchanged Mf^ ~ 984 MeV, Cfg-KTr increases from 1.35 to 2.35 GeV^^and the ^-^ ^ increases from 0.9 to 2. 
In Fig. [3] the dotted line corresponds to the case Cf^KK = 12 GeV~^ of Table [Til while the data above 0.8 GeV 
are well described, an overestimation of data point between 0.6 and 0.8 GeV is clearly visible. Thus, in the present 
analysis a destructive interference with a meson should occur: this fact can represent a constraint on models of 
light scalar mesons. Interestingly, a destructive interference of /o(600) and /o(980) is also the outcome of Ref. [2l] |. 
where the /o(980) channel overestimates the data below ^ 700 MeV. Nevertheless a more refined analysis including 
the interference with the /o(600) meson and using the new data of KLOE [34j -for which at present no tables with 
branching ratios have been presented- will be compulsory and represents an outlook of the present work. At this stage 
a preliminary comparison with the new data of KLOE can be done only by considering the integrated branching ratio 
given m Ref. H: BR{(t> 7r°7r°7) = 1.07 x lO"". Our value runs from 1.06 — 1.10 x 10 depending on the choice 
for CfgKK, and therefore is consistent with the latest experimental results. 

We now compare and discuss our amplitudes extracted from the KLOE and SND data with the amplitudes extracted 
from the experimental analyses of ^] for the /o meson via j/ip decay at BES. Let us stress that a comparison by 
using the decay widths would be less reliable, since they depend on the adopted way to evaluate them. The amplitudes 
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TABLE II: Fit of the (ji ^tx^-k'^ line shape data from x > 800 MeV. 



CfoKK 


(MeV) 


C/oTTTT 




2 

X 

d.o.f. 


(GeV) 


(GeV) 





984.2 ± 3.2 


1.35 ± 0.04 


0.263 ± 0.007 


0.9 





1.39 ± 0.04 


2 


984.1 ±3.0 


1.38 ± 0.04 


0.275 ± 0.007 


0.9 


0.68 ±0.01 


1.42 ±0.04 


4 


984.0 ± 2.8 


1.48 ±0.05 


0.309 ± 0.009 


1.0 


1.36 ±0.01 


1.52 ±0.05 


8 


983.6 ± 2.4 


1.81 ±0.07 


0.419 ± 0.013 


1.5 


2.72 ±0.03 


1.86 ±0.07 


12 


983.4 ±2.4 


2.25 ±0.09 


0.558 ±0.019 


1.9 


4.07 ±0.04 


2.32 ±0.09 


14 


983.3 ± 2.3 


2.50 ±0.11 


0.632 ±0.02 


2.0 


4.75 ±0.04 


2.57 ±0.11 



extracted in Refs. [13, are : 

A-f„^^ = 2.88 ±0.22 GeV, Af^KK = 5.91 ± 0.77 GeV. (27) 

Our last entries in Tables HI and [TTl are in qualitative agreement with Eq. (|27p . A more quantitative check is possible: 
we deduce from Eq. (|27[) the couplings Cfg^rTr and CfgXK (at a fixed value of the mass of the meson) and use them 
in the fit of the KLOE and SND data leaving the coupling c^^jg as the only free parameter. The fit turns out to be 

acceptable with j^-^^ ~ 2.4, the extracted values of the couplings are c/qtttt = 2.79 GeV~^ and cj^kk — 17.37 GeV~^ 
for M/q = 984 MeV. The fact that the line shape of (/) ^ ^tt^tt^ can be well described by using Eq. ((27)) as an input 
shows that the BES results of Ref. [1^ are compatible with our analysis. 

B. The (f) -» -y-K^r) fit 

We proceed as in the /o case by first setting CagKK = 0: the fitted curve is shown by the solid line in Fig.[4l We then 
increase CaoKK up to 14 GeV~^: the fit turns out to be satisfying in the whole range, see Table ITTTl where the mass, 
the couplings and the on-shell amplitudes Aag-Kn = ^ao7r?)(^ao), Aa^KK = AagKKiMag) defined in Eq. (PT|l - (|^ . are 
reported. In Fig. [4] we also show the fitted curve in the case CagKK = 12 GeV~^ (dashed line). The narrow peak due 
to threshold effects is remarkable, but the data are not precise enough to determine its existence. 

2 

A comment about the behavior of the ^'^ ^ as CagKK varies is in order: the fact that the fit is quite good for 
all the values of CaoKK is a signal of a strong correlation between the parameters. As in the /o case CagKK only 
appears in the denominator of the propagator, thus the dependence of the line shape on this parameter is weak and 
its determination hard. While the value of the mass is almost constant, CagTrri and c^^ao change significantly. This 
means that, while the fitted curve is in agreement with the data, the determination of the parameters is not reliable. 
This is due to the large error bars of the data, which in turn produce also large errors associated with each of the 
parameters of the fit. Again a new analysis with the new data of KLOE would be extremely important. As before 
a comparison with the new results of KLOE on the integrated branching ratio is useful. In Ref. [13] it is found 
that BR{(f) 7777^7) = 6.92 — 7.19 x 10~^ (respectively for the chains rj ^ and rj 7r+7r~7r°). Our value runs 
from 7.37 — 7.73 x 10~^, depending on the choice for CagKK, and therefore is in reasonable agreement with the latest 
experimental results. 

We now compare our results with the amplitudes extracted from the experimental analysis of tlie Crystal Barrel 
data of Ref. 

Aag^r, = 3.33 ± 0.15 GeV, AagKK = 3.59 ± 0.44 GeV, (28) 

which are not far from the last three entries of Table lllll As done previously, a more quantitative test consists in 

deducing the couplings CagTiri and CagKK from Eq. (j28p and, at a fixed value of the mass, perform the fit to SND and 

2 

KLOE data with only c^^ap as free parameter. The fit turns out to reproduce correctly the data with ^-^ ^ ~ 1.2, 

the extracted values of the coupHngs are Cag-Ki^ = 10.17 GeV^^ and CagKK ~ 9.79 GeV^^ for Mag = 1004 MeV. Also 
in this case we obtain that our analysis is compatible with experimental results of a different source. 
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C. Comparison with the fit without derivatives 

As we have seen a structureless approach using derivative couplings describes well the KLOE and SND data of the 
(j) radiative decays. This is manly due to the fact that within the derivative approach one has an extra dependence of 
the form (x^ — m\ — m^)^ in the numerator of the line shape, which emphasizes the peak of the distribution at large 
values of the invariant mass. Since this term is not present in the non-derivative case the description of the peaked 
line shapes more difficult. As a comparison we also discuss in the following the results obtained when non-derivative 
interactions are employed. One could still obtain the peaks with non-derivative couplings by using very small decay 
amplitudes, which however are not in agreement with the large amplitudes extracted from Rcfs. 24, 25], reported in 
Eqs. (EZI) and (Hg). 

More specifically, the inconsistency of the use of non-derivative couplings is evident in the case of the /o meson: 
when setting Cf„KK = c/qth- = and allowing for nonzero df„KK and d/(,,r7r (see discussion in the Introduction) high 

2 

J >12 are obtained for fits above 800 MeV. 
When setting in the ao channel CaoKK = Ca„TTri — and allowing for nonzero daoKK and daoirr] fits with still 

2 

acceptable values of j ~ 2.4 are obtained, but the resulting amplitude Aaa-nr^ ~ 1.37 is smaller than all the cases 
corresponding to the derivative coupling (see Table [1111) and therefore the comparison with the amplitudes of Eq. 
is problematic. Reversing the argument, we may deduce the values of dag-Kri and da^KK from Eq. (p8)) (again fixing 
the mass of the meson) and use them to fit the KLOE and SND data with one free coupling c^-yag . The obtained 

2 

line shape describes very badly the data with ^ ^ 10, much larger than all the cases corresponding to derivative 
coupling. Thus, incompatibility with Eq. (|28p is manifest. 

All these results indicate that the use of non-derivative coupling is disfavored by present experimental informations. 



IV. CONCLUSIONS 



In this work the role of derivative interactions of the mesons /o(980) and ao(980) with the pseudoscalar Goldstone 
bosons, which are a basic consequence of spontaneous chiral symmetry breaking as shown explicitly in chiral pertur- 
bation theory, has been investigated within a nonlocal approach in relation to the radiative decays ip ~* jn'^Tr'^ and 
(j) 'yn'^ri. After developing the method at one-loop level it has been shown that a satisfactory description of the 
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TABLE III: Fit of the (j> -> 77r°r? line shape. 



CaoKK 


Moo 
(MeV) 


CaQ 77 77 


(Gel/-i) 


d.o.f. 


AagKK 

(GeV) 


(GeV) 





1005 ± 18 


4.68 ± 0.90 


0.263 ± 0.046 


0.88 





1.54 ± 0.3 


2 


1005 ± 18 


4.71 ±0.94 


0.278 ± 0.050 


0.88 


0.74 ± 0.05 


1.55 ±0.32 


4 


1004 ± 18 


4.80 ± 1.07 


0.321 ±0.06 


0.88 


1.47 ±0.10 


1.57 ±0.36 


8 


1003 ± 18 


5.10 ± 1.62 


0.486 ±0.13 


0.88 


2.92 ±0.20 


1.67 ±0.53 


12 


1002 ± 19 


5.52 ± 2.49 


0.731 ±0.28 


0.92 


4.36 ±0.31 


1.80 ±0.82 


14 


1001 ± 20 


5.77 ±3.01 


0.873 ± 0.400 


0.92 


5.09 ±0.37 


1.88 ± 1.00 



line shapes of the reactions cf) — > jtt'^it'^ and (f> ~^ jTr'^'i] is obtained: the property of derivatives interactions seems 
tailor-made to describe the well-marked peaks measured experimentally. 

More specifically, in the (f) —> 77r''7r'^ case we fitted our theoretical curves, which make use of derivative interactions 
only and involve a virtual /o meson, to data above 600 and 800 MeV (Fig. [3] and Tables U and |TJ) : in the latter 
case the influence of the /o(600) resonance is less relevant. It has been also stressed that the determination of the 
coupling to kaons from the line shape only is hard because the latter enters only in the propagator of the scalar meson. 
Interestingly, fits involving a large /o KK coupling, as found in [l^l, are acceptable (x^/ d.o.f. < 2) when data- 
points above 800 MeV are considered. Notice that in this case data below 800 MeV are overestimated as the dotted 
line in Fig. [3] shows. A destructive interference in the /o(600) channel should be invoked in order to get agreement 
with data, see also Ref. [2l|. While leaving a more detailed study as an outlook, this fact can represent a significant 
constraint on models based on a particular interpretation of the scalar mesons, such as the tetraquark assignment. 
We also verified that the usage of non-derivative couplings generates a bad agreement with the experimental results 
{x^ /d.o.f . > 10): the line shape of the /o meson cannot be reproduced for any choice of the parameters. 

In the ^tt'^t] case, occurring via a virtual oq meson, a good agreement is found (Fig. |4]and Table [TTT|) . However, 
even sizable variations of the parameters do not make the agreement worse: for this reason a determination of the 
parameters is hard. As above, the determination of the coupling to kaons is problematic. However, when using the 
amplitudes extracted from Ref. f23| to deduce the couplings with pseudoscalar mesons a good description of the line 
shape is obtained: this fact shows compatibility of the study of Crystal Barrel data in Ref. [2^ with KLOE-SND data 
when derivative interactions are used (notice also the pronounced peak close to threshold of the dashed line in Fig. 
m whose appearance is due to a large coupling to kaons). On the contrary, the same procedure in the non-derivative 
case shows that the line shape of (/> ^ 777^77 is badly reproduced. 

Summarizing, the results of the present study point out that derivative-type interactions (denoted as mechanism 
a.l in the Introduction) with pseudoscalar mesons can play an important role for the study of the scalar states ao(980) 
and /o(980). The next step of the present work is the inclusion of non-derivative interactions (explicit chiral symmetry 
breaking terms, mechanism a. 2) besides the derivative ones. The next-to-next (more ambitious) step is the inclusion 
at the same time and in the same theoretical framework, together with the mechanisms a.l (derivative couplings) 
and a. 2 (non-derivative couplings), also of the correspondent kaon-loop driven contributions (mechanisms b.l and b.2 
respectively). To achieve these goals, which can also set up which mechanism (if any) is dominant in radiative (j) 
decays, further work is required. The aim is a definitive determination of the relevant amplitudes, which represent a 
necessary tool to test models and thus to understand the nature of the light scalar mesons. To this end the use of other 
experimental informations from j /ij} decay and Crystal Barrel data is necessary. Also, a comparison with the newest 
KLOE data [13] and the study of radiative decays of scalar resonances [35| represents an interesting development. 
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